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Abstract. Magnetic flux tubes reaching from the solar convective zone into the 
chromosphere have to pass through the relatively cool, and therefore non-ideal 
(i.e. resistive) photospheric region enclosed between the highly ideal sub-photospheric 
and chromospheric plasma. It is shown that stationary MHD equilibria of magnetic 
flux tubes which pass through this region require an inflow of photospheric material 
into the flux tube and a deviation from iso-rotation along the tube axis. This means 
that there is a difference in angular velocity of the plasma flow inside the tube below 
and above the non-ideal region. Both effects increase with decreasing cross section 
of the tube. Although for characteristic parameters of thick flux tubes the effect 
is negligible, a scaling law indicates its importance for small-scale structures. The 
relevance of this "inflow effect" for the expansion of flux tubes above the photosphere 
is discussed. 
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1. Introduction 

The interaction of solar flux tubes with the surrounding plasma is usu- 
ally treated in the framework of ideal MHD (i.e. with zero resistivity), in 
which no exchange of plasma between the flux tube and its environment 
is possible. While this approach appears to be well suited for both the 
convection zone and the upper chromosphere (where the degree of ion- 
isation is sufficiently high), it becomes doubtful for the relatively cold 
and therefore almost neutral photospheric plasma (see Figure 4) . It this 
resistive layer, deviations from the rigid coupling between fluid and field 
must be anticipated. This could have important consequences for the 
widely used conception of flux tubes being wound up by photospheric 
motions. Also the strict separation of plasma within the flux tube from 
its environment as required in ideal MHD might break down in this 
resistive layer. The purpose of this work is to compute this deviation 
from ideal MHD in a self-consistent manner. 

First we consider a stationary magnetic flux tube with both ends an- 
chored in the convective zone (see Figure 2a). The flux tube can be 
thought of as consisting of a set of nested tubes which are flux surfaces 
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Figure 1. Orientation of unit vectors e p and e t with respect to the tube surface and 
the lower photospheric boundary. 



for the magnetic field. Assigning to each of these surfaces the magnetic 
flux it encloses defines a function Vl/(r), B(r) • = 0, which is zero 
on the tube axis and monotonously increases outwards. (We assume 
that there are no field reversals within the flux tube.) In a stationary 
situation, any ideal MHD flow has to preserve these flux surfaces, and 
the plasma velocity has to be tangential to the surfaces of constant 
v • = 0. The lower boundary of the domain under consideration 
(given by the lower boundary of the photosphere) is a surface which 
intersects the flux tube twice. Any plasma motion imposed on the 
boundary at one footpoint implies a corresponding motion at the other 
footpoint. The exact relation between these motions is derived from 
"ideal" Ohm's law (i.e. as it is known from ideal MHD) 



At the lower boundary, the plasma velocity and the magnetic field may 
be decomposed into their poloidal and toroidal components: 



The toroidal components are directed along the intersection of the 
boundary surface with the ^-surfaces of the flux tube. Their orientation 
can be defined by requiring that the toroidal unit vector e t have a 
positive orientation with respect to the magnetic field vector on the 
tube axis. The poloidal components are also tangential to the ^ surfaces 
but perpendicular to e t , with a unit vector e p := x e t / ||V\I/|| (see 
Figure 1). Assuming v p = and using B p = b p V^> x e t and v t = ut e t , 
Equation (1) yields 



-V$ + v x B = . 



(1) 



v = v t + v p and B = B t + B p . 



(2) 



v t x B p = v t b p V* 
1 d$ 



(3) 



(4) 
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since (3) implies $ = <l>(^). This equation shows that a given v t distri- 
bution at one end of the flux tube determines d&/d*k, a function which 
depends only on ^ and is thus constant along the flux surfaces, thereby 
inducing a corresponding v t distribution at the other end. For a flux 
tube which is perpendicular to the boundary and which has circular 
flux surfaces ^ = VE^r) (where r is the distance from the tube axis), 
the poloidal component of B is given by B p = l/(27rr) V\&(r) x e t and 
hence b p = l/{2nr). In this case ft and the angular velocity oj := v t /r 
are functions of r only, i.e. they are constant on each flux surface This 
is simply Ferraro's law of iso-rotation (Moffat, 1978). In the general case 
vt is not constant on flux surfaces, but an integration along e t yields 
the circulation time 



which only depends on the flux surface. This quantity (or the corre- 
sponding angular velocity = 2ir/T) explicitly shows the coupling of 
the toroidal velocity field between both ends of the flux tube. 
While the preceding results were based on the idealness of the plasma, 
we will now investigate the effect of a non-ideal region the flux tube 
has to pass. This non-ideal region is given by the comparatively cold 
photosphere. Here a possible slippage effect due to the non-ideal pho- 
tospheric region would result in a deviation of v p from (4) . Also in the 
case of incompatible poloidal velocities on both footpoints the onset of 
slippage will keep the resulting twist of the flux tube finite, as opposed 
to the infinite "winding-up" of field lines expected for ideal MHD. 



To study the effect of a resistive layer on the flux tube, it is sufficient to 
consider only one half of the tube and concentrate on the photospheric 
region close to the footpoint, as shown in Figure 2a. For simplicity, 
we will restrict ourselves to stationary, axisymmetric solutions. The 
ensuing calculations will use cylindrical coordinates [r, (f), z], with unit 
vectors [e r ,e^,e z ]. The (z = 0)-plane is given by the photosphere's 
lower boundary, while the z-axis coincides with the tube axis and is 
pointing away from the Sun. The problem's axial symmetry is now 
conveniently incorporated by setting d<j, = 0. With dt = 0, the set 
of MHD equations to be solved for the mass flow velocity v and the 
fields B and E = — consists of the momentum balance (6), a resis- 
tive Ohm's law (7), the equation of continuity (8), and the remaining 





2. The model 
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Maxwell equations (9, 10): 



= -VP + j xB + pg 

77 j = -V$ + v x B 

= V • (p v) 

fj, j = V x B 

= VB 



(6) 
(7) 
(8) 
(9) 
(10) 



As usual, p and 77 denote the plasma's mass density and resistivity, 
respectively. The inertia term p(v ■ V)v is omitted from (6) since its 
ratio to the induction term j x B is of order 0[(v/va) 2 }, where 



is the Alfven velocity. Adopting B ps 0.1 T and p ps 10 -6 kg m~ 3 as 
characteristic values for our photospheric flux tube yields v\ « 90 km 
s _1 , which is large compared to the magnitude of observed photospheric 
plasma motions of v ^, s ps 5 km s _1 . Section 6.3 gives an a posteriori 
verification of this conjecture. 



Figure 2. Left: Flux tube emerging from the Sun's photosphere. The white box 
marks the section on which our computations focus. Right: Our flux tube model. 
The shaded area indicates the poloidal plane in which B p , vi and e^j. are located. 



3.1. Derivation of the Inflow Equation 

An important difference between the ideal and the non-ideal case is 
the exchange of plasma between the flux tube and its environment, a 
process which is impossible in ideal MHD. The plasma flow across the 




(11) 




Figure 2 a. 



Figure 2 b. 



3. Resistive Inflow 
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flux surfaces of the magnetic field can be derived from Ohm's law (7) 
alone 

-V$ + v x B = fj V x B (12) 

with (9) inserted and fj := rj/fi substituted. Again we decompose v and 
B similar to (2) into their toroidal and poloidal components, where the 
toroidal component is directed along and the poloidal plane is the 
r — z-plane. This yields, after insertion into (12), 

v p x Bp = fj V x B p (13) 

as the poloidal component of Equation (12). Now let 

ey := -j- Bp and e± := e t x ey (14) 

-Dp 

be two orthonormal vector fields parallel and perpendicular to B p . Then 
the crossproduct of Equation (12) with B p , together with v p = v±ej_ + 
vyey yields: 

Bp x (v p x Bp) = ?} Bp x (V x B p ) 
O v± e ± = fj Bp x (VB p x ey + B p X7 x ey) 

v± e ± = fj (e ± -V(ln| J B p |)-(Vxe||)-e t )e ± . (15) 

3.2. Discussion of Inflow Properties 

From the "inflow equation" (15), the following flow properties are 
evident. First, the flow magnitude is proportional to fj and thus, as 
expected, vanishes as soon as ideality is restored. Second, both mag- 
nitude and direction of B p play no role for v_|_ or, in other words, a 
substitution B p — > ±a B p leaves vj_ unchanged for any constant a. 
(Note that a substitution B p — > — B p changes the direction of both en 
and ej_, thereby preserving the direction of Vj_.) 

Moreover, since for a flux tube B p generically decreases outwards, there 
will generally be an inflow of matter into the tube throughout the 
entire region where rj ^ due to the first term on the rhs of Equation 
(15). The contribution of the second term will be negligible in generic 
cases for the following reason. Both terms define a characteristic length 
scale. For the first term, this is the scale R on which the poloidal field 
decreases markedly. It can be used for defining the radius of the flux 
tube as well. The second term defines a typical curvature radius R c 
of the poloidal field lines. If R c is of the order of R, the flux tube is 
strongly distorted, i.e. the change of its cross section is of the same 
size as the cross section itself. A closer analysis shows that the field 
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lines of B p have to be bent strongly inwards for the second term 
to contribute to an outward directed flow and to dominate over the 
first term. (For an instructive example see the Appendix.) However, 
observational evidence suggests that a flux tube's cross section either 
stays more or less constant (Klimchuk (2000),Watko and Klimchuk 
(2000)) or increases monotonously with height, as in sunspots. Notice- 
able amounts of inward curvature are produced in neither of the two 
cases, and the second term of (15) can thus be ignored without much 
loss of generality. (Note that even if such cases should occur, the notion 
of a tube-shaped configuration requires that strong inward curvature 
of field lines at one tube part be balanced by a suitably strong outward 
curvature at some other part. Consequently, the weaker the inflow gets 
at one point, the stronger it gets at some other point, as can clearly be 
seen in Figure 8 of the Appendix. Although the fact that the net value 
of this mutual cancelling of flow depends on the global density structure 
makes a precise quantitative treatment of the most general case more 
difficult, it seems reasonable to assume that even then the net inflow 
will be diminished only moderately by strong poloidal curvature.) 
In the case of straight flux tubes, the approximation of small V x ey 
becomes exact and leads to a scaling of the inflow velocity 

||v_l|| oc l/R , (16) 

which means that the inflow is more violent for thinner tubes. For 
instance, comparing cylindrical flux tubes with the same B p profile 
but different characteristic radii R 



B p (r) = B b z (r/R) e z (17) 

we find 

v ± = f, V(ln \B P \) = f) R' 1 8 x (\n \b z (x)\) (18) 

where the dimensionless radial coordinate x := r/R has been intro- 
duced. We may thus conclude that the total mass inflow through a 
cylindrical surface of radius R occurring within the resistive layer, 



:= J(pv ± )-da = J(pv±) 2nRdz (19) 
(d x ln|M*)l) J p(z) v(z) dz , (20) 



M 

= 2tt 



is scale-independent with respect to R under the assumption of a hor- 
izontally stratified atmosphere (p = p(z) and fj = fj(z)), i.e. tubes of 
various radii but with the same b z profile will transport the same mass 
rate, regardless of their strength. (Note that the momentum equation 
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(6) was not used to derive the preceeding results, which therefore are 
not limited to flow fields satisfying v <C va-) 

4. Solving for the Complete Flow Field 

From now on, it will be assumed for simplicity that the tube is of 
cylindrical shape, i.e. it does not "fan out" . This simplifying assumption 
is justified by the fact that the main effects of the resistive layer, namely 
the flow of plasma into the flux tube and the decoupling of toroidal 
velocities above and below this layer, are both already present in this 
simplified geometry. Our assumption then translates to B r (r, z) = 0, 
so that the solenoidality condition (10) now reads 

= V B = d z B z (r,z) . (21) 

Now consider the momentum equation (6) in its [r, (f>, z] components: 

/0\ fd r P(r,z)\ ( \ 

=- +j xB + p(r,z) (22) 

\0J \d z P(r,z)J \-gJ 

From the ^-component we derive d z B^{r, z) = and hence B = B(r), 
such that we can integrate the pressure from the r-component of (22) 

P(r,z)= Pl (r)+p 2 (z) , (23) 

which in turn leads to d r p(r,z) = due to the z-component of (22). 
Assuming that P, p and T are linked by an arbitrary equation of 
state P = P[p,T], a horizontally stratified temperature T{z) yields 
d r P(r, z) = and hence 

j x B = . (24) 

Therefore a cylindrical flux tube has to be force- free if the temperature 
of the plasma is horizontally stratified. 

Although the requirement of strictly horizontal isotherms T = T{z) 
is clearly not fulfilled for very thick tubes (as readily seen by the 
reduced intensity observed in sunspots), we choose to adhere to this 
assumption not only for the benefits of a markedly simplified analyt- 
ical treatment but also for physical reasons. Setting d r T = would 
be justified provided that sufficiently strong horizontal heat transport 
was present. However, under the assumption of an ideal plasma, the 
(radiation-dominated) influx of heat does not suffice to heat the interior 
of an embedded flux tube to the ambient temperature level because 
convective energy transport is inhibited by the tube's strong magnetic 
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Figure 3. Selected field lines of = (x/{\ + x 2 ), 1/(1 + a; 2 )) , which satisfies 

(25) and corresponds to G(x) — [2(1 + a; 2 )] -1 . 

field, as was first realised by Biermann (1941). But within the non- 
ideal zone, convection across magnetic surfaces is well permitted (or 
even enforced, see Section 3), such that the radial exchange of heat 
will be amplified significantly. The inflow effect thus reduces the radial 
temperature gradient and may possibly lead to thermal structures with 
d r T <C d z T, in which case our assumption would be clearly justified. 
(This condition should be easily fulfilled for thin tubes, while thick 
tubes will hardly be affected by this reasoning since it takes too long 
to exchange noticeable fractions of their mass contents via the inflow 
effect, see Section 6.2.) Moreover, the presence of neutral gas in this 
region allows for additional convective energy transport unimpeded by 
magnetic fields. 

In our case, the only non-trivial component of (24) is the r-component, 
which reduces to 

d x {b 2 z (x) + bl(x)) +2 6j(x)/x = 0, (25) 

where again x = r / R, b<f,(x) := B^xR) / Bq, and b z (x) := B z (xR) / Bq 
were used. This equation is solved by 

b 2 ( x ) = _ x d x G(x) and b 2 z {x) = (1/x) d x [x 2 G(x)} (26) 

for any function G(x) satisfying d x G(x) < and d x [x 2 G(x)] > 
(Schliiter, 1957). Figure 3 shows a typical solution to (25). 
Given a force-free magnetic field, the solution for the flow velocity v 
is to be determined from Ohm's law and the equation of continuity. 
This also requires to fix boundary conditions for the velocity on either 
the upper or the lower boundary of the domain. Here the linearity of 
the two equations with respect to v is very useful because any solution 
can be seen as a superposition of a solution v;d of the ideal Ohm's law 
(1) and the continuity equation, and a particular solution v res of the 
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resistive Ohm's law (7) and the continuity equation: 



v = v id + v res . 



(27) 



Since we are interested in the deviation from iso-rotation due to the 
resistive photospheric layer, we are free to determine v rcs for a certain 
choice of boundary conditions, and the solution for any other boundary 
condition is then given by adding a corresponding ideal solution. The 
most simple boundary condition is setting v = on the upper boundary 
z = z up , such that the velocity on the lower boundary exactly equals the 
difference of the toroidal velocities above and below the photosphere, 
i.e. the deviation from iso-rotation. In this case we have (the index 
"res" on the solution is suppressed in the following): 



v r (x, z) = —R 1 P\{x) rj(z) 
V(t> {x,z) = ±R~ 2 [/3 2 (x) I 2 {z 
v z (x,z) = +R- 2 f3 4 (x) I 2 {z) 



+ h{z) 



(28) 



where the Pi..a(x) are given by 
b'Jx) 



ft 0*0 : = 
h(x) := 



+ 



b z (x) 

b<f,(x) 
b z (x) 

bcfyjx) 
b z (x) 

b'4x) 
b z (x) 

b" z (x) 



b" z (x) + 1 b' z (x) 



x b z {x) 
b'Jx) 



+ 



b z (x) 

b" z (x) 
b z {x) 

K(x) 

b z{x) V 6 (x 
1 b' z (x) 



+ 



+ 



(29) 
(30) 

(31) 



b z (x) xb z (x) 
and the h^z) are defined as 



K(x) 
b z (x) 



h(z) :- 



r/(C) dC and I 2 (z) 



ri(0 



p(0 



(32) 



(33) 



The sign of v ( f,(x, z) in (28) is opposite to the sign of b^(x), which is not 
fixed by (25) and may be chosen arbitrarily. (Figure 3 has sgn b^ = +1.) 
To proceed further, one could either prescribe a vortex at z = and 
use (25) and (28) to compute the magnetic field components, or insert 
into (28) a typical solution of (25). The first alternative would use a 
(rather long and messy) first order ODE, while in the latter case the 
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flow field could be read off directly from (28). Therefore, this avenue is 
chosen here. 



5. Realistic Input Parameters 

For a quantitative evaluation, prescription of density and resistivity 
profiles p(z) and rj(z) is required. We use the data provided by the solar 
atmosphere model "C" of Vernazza, Avrett, and Loeser (1981) (here- 
after VAL, see Figure 4) along with the conductivity calculations by 
Kubat and Karlicky (1985) based on the VAL model. Since this model 
neglects magnetic forces, we continue to assume isotropic resistivity for 
both simplicity and consistency, that is, we take r/ = (en) from Table 
III in (Kubat and Karlicky, 1985) (see Figure 5). The function 

-l 

(34) 

with [770 = 0.058 O m, z m = 360 km, L v = 330 km] also depicted there 
will be used to model the photosphere's actual resistivity; its density 
is approximated by 

p(z) = po exp(-z/L p ) (35) 

with po = 3 • 10 -4 kg m -3 and L p = 120 km, which is in sufficient 
agreement with the VAL model data for p(z). 



8000 - 



5000 - 




4000 500 1000 1500 2000 

z [km] 

Figure 4- Temperature variation with height according to the VAL model. The 
shaded area marks the resistive region, in which the ionisation ratio drops below 
10 -5 . The minimum of T(z) at z ~ 500 km is clearly discernible. 



rj(z) 



1 + 



z Z B 
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Figure 5. Computed resistivity profile (crosses) vs. analytic model function. The 
relatively poor agreement between relation (34) and the data at z < is of minor 
importance since these layers are not explicitly considered here. 

In the ensuing quantifications, we will specialise to the B field of Fig- 
ure 3 as a "flux tube prototype". This seems justified since tentative 
computations using other fields have yielded only very small deviations. 
Additionally, z up — > oo is used since above the non-ideal region the 
contribution to (33) becomes negligible. 
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6. Quantitative Flow Evaluation 

6.1. The Scaling Law 

According to (28), there must be a tube radius R tr . such that 

— ii ii / R 1 '• R S> i?tr. toa \ 
V = II V H * { iT 2 : R«R tl . ^ 

Inserting our parameters found in the preceeding section, we find R tr . ~ 
5000 km. Since at this radius v will have dropped below 1 m s _1 , we 
may safely regard 

||v||ocir 2 (37) 
as the relevant scaling law for small scale flux tubes. 

6.2. The Blowup Timescale 

The knowledge of absolute photospheric density and resistivity allows 
us to quantify the total mass inflow (19) associated with a cylindrical 
tube as M « 1.9 • 10 8 kg s~ 1 , which, when compared to the total mass 

M« := J „ dV = .*> J" p(z) d 2 « 4.9 ■ 10" kg (j^) 2 

(38) 

of the plasma contained inside the tube, defines a typical timescale 

rz U p 

Mtot R 2 L P{z)dZ / R \ 2 

blowup := = 70 h — — (39) 



M 2 P up ( \ *• ( \ j ' '"" V!00 km 
Jo 

at which the tube exchanges a noticeable fraction of its contents. 
6.3. Condition for sub-Alfvenic Flows 

Since the flow magnitude scales oc R~ 2 , the requirement v <C v\ is 
actually a limitation on the radius of the flux tube. According to Figure 
6, this may be quantified as R Z 10 km, which is well below the re- 
solvable scale achieved by present (and near-future) solar observations. 
Discarding the inertia term in (6) was indeed justified. 
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6.4. Field Line Slippage 

The toroidal flow depicted in Figure 7 shows a striking deviation from 
the flow expected in the ideal case rj = (in which Ferrano's theorem of 
iso-rotation forces all field lines to rotate with constant angular velocity 
Q(x,z) := v$(x, z)/(xR), such that d z Q, = 0). Note also that far below 
the photosphere (z < 0), the lines of constant tend to become 
vertical, implying that iso-rotation is recovered as r\ again tends to 
zero. 

Since we have l2(0)/Ji(0) ~ 0.1 <C 1, the profile of the footpoint vortex 
(i.e. the cut along z = 0) can be approximated by 

M x >°) _ ( R \ ( 40 ) 



130 m s- 1 V 100 km / 1 + x 2 
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With Q(x, z up ) = 0, the total difference in angular velocity below and 
above the non-ideal layer is given by 



For observable tube sizes (R <; 100 km) this would require a velocity 
resolution close to (Av ) reso := R AQ, ~ 30 m s _1 . Although this limit is 
not quite reached by current imaging techniques, the further improve- 
ment in image resolution may soon render observational verification 
feasible. 



7. Implications for the Tube's Global Evolution 

Since, according to our results, plasma has to flow into the tube from 
both ends and cannot leave the tube outside the non-ideal zone, the 
question arises as to where the inflowing matter goes. The possibilities 
are a) a steady increase of the tube's volume (tube gets "inflated"), b) 
tube is static and downflow into the convection zone occurs or c) the 
inflowing plasma recombines within the flux tube and leaves the tube 
in the form of neutral gas. (Of course, in reality various combinations 
of a) to c) are conceivable.) In the present model, the direction of 
vertical flow inside the tube is determined by the boundary condition 
at z = z up (or any other height), such that up- or downflows of arbitrary 
magnitude can be achieved by choosing a correspondingly large (pos- 
sibly negative- valued) profile for x i— >■ v z (x, z up ). However, we have no 
reason to favour any specific profile, and thus our present, rather simple 
model cannot provide a definitive answer here. (Note that in Section 4, 
v| 2up =0 (implying downflow at z = 0) was merely chosen to simplify 
the calculation of field line slippage. It was not supposed to indicate a 
preference for downflows in any way.) To shed light on this important 
issue, the aforementioned possibilities a) to c) suggest two avenues for 
an extension of our model. First, if plasma was flowing up the tube, 
thereby forcing it to expand in length and/or cross section, the tube's 
field lines would be stretched, and their tension increased. Eventually, 
the growing contribution from the j x B force might become strong 
enough to balance the gas pressure, causing the upflow to cease. To see 
whether such a final equilibrium state exists, and if so, what the tube 
parameters in such a state are, one would need to abandon cylindrical 
symmetry and model the full arch-shaped tube such that all effects of 
field line curvature could properly be accounted for. Unfortunately, the 



Afi(x) : 




< 



AO(0) « 3.1 • 10~ 4 s- 1 [12/(100 km)]" 3 . (42) 
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corresponding set of equations could turn out to be very difficult to 
solve analytically, and thus the feasibility of this approach is unclear 
at the moment. 

Second, evaluating the significance of possibility c) would obviously re- 
quire the introduction of radial gradients of ionisation and temperature. 
In such a model, one-dimensional reference atmospheres such as VAL 
can no longer be used to prescribe atmospheric parameters (except at 
large distances from the tube axis), and self-consistent modelling of 
density and temperature becomes mandatory. Again, it seems doubtful 
whether analytic solutions can be obtained at a reasonable expenditure. 
Still, in both cases a recourse to numerical investigations of the de- 
scribed settings remains a vital option and may help to clarify the 
role of the inflow effect with respect to the tube's global temporal 
evolution. (A discussion of observational evidence for downflow is given 
by Frutiger and Solanki (1998), but whether these observations can be 
applied to the photospheric region remains unclear since they refer 
to velocities measured at coronal or transition region temperatures. 
Generally speaking, the very existence of pronounced vertical flows 
inside photospheric flux tubes still seems to be a controversial issue 
among the observing community.) 



8. Summary 

Our analytic investigation of stationary MHD equilibria of magnetic 
flux tubes has shown that Ohm's law enforces an inflow of fluid to- 
wards loci of higher field strength, which depends neither on the tube's 
cross section, nor on the strength and direction of its B field. Being 
proportional to rj, this inflow occurs wherever the tube penetrates the 
cool photospheric layer, in particular at the tube's footpoints. 
It was shown that a static flux tube of cylindrical shape has to be 
force-free if the ambient plasma temperature is horizontally stratified, 
a result which holds for arbitrary values of plasma beta. The intro- 
duction of a resistive layer allows for stationary MHD solutions with 
finite field twist and a difference in the rotational velocity above and 
below this resistive layer. This constitutes a marked deviation from the 
iso-rotational behaviour known from ideal MHD and limits the winding- 
up of the flux tube's field lines if incompatible rotational velocities are 
imposed on the tube's footpoints. Although according to the scaling 
law for the plasma flows these effects either too small or too slow to be 
detected by present solar observations (i.e. the effect either requires too 
small structures or produces velocities below the detection threshold) , a 
future improvement of observational resolution may soon show whether 
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the described effects can be distinguished from the convective motions 
of the ambient plasma. 



Appendix 

To clarify some aspects concerning the direction of radial flow, consider 
the magnetic field 



B 



sp 



Bq coshz 
1 + (r coshz) 



— r sinhz, 0, coshz 



(43) 



where the r and z coordinates are now dimensionless for simplicity 
of the argument. The tube radii have a (cosh z)' 1 profile and the 
field decays as ||B sp || z= o oc 1/(1 + r 2 ). Figure 8 shows a vector plot 
of the corresponding perpendicular flow component vj_. Beyond the 
dotted line, the field curvature gets so strong that the flow direction 
is indeed reversed, leading to an outflow of plasma. Although the ex- 
istence of photospheric fields like (43) cannot be ruled out completely, 
the rather low astrophysical significance of such configurations (as dis- 
cussed in Section 3.2) is further diminished by the fact that B sp has 
V x ( j x B) sp / and therefore does not describe a flux tube in the 
hydrodynamical sense. 




Figure 8. Vector field plot of vi and selected field lines of B (solid) in the poloidal 
plane. The flow direction reverses at the dashed line. Note how high curvature 
increases both inflow (due to the "inflow effect" ) and outflow (due to the requirement 
that the tube be void of field line reversals). 
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